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ABSTRACT 
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I ■ We present "On The Fly" maps of the CO(l-O) and C0(2-l) emission covering a 

' 10' X 10' region of the NGC 6946. Using our CO maps and archival VLA HI observa- 

. tions we create a total gas surface density map, Soas, for NGC 6946. The predominantly 

molecular inner gas disk transitions smoothly into an atomic outer gas disk, with equiv- 
alent atomic and molecular gas surface densities at R = 3.5' (6 kpc). We estimate that 
\ the total H2 mass is 3 x 10^ Mq, roughly 1/3 of the interstellar hydrogen gas mass, 

' and about 2% of the dynamical mass of the galaxy at our assumed distance of 6 Mpc. 

^ ■ The value of the CO(2-1)/CO(1-0) line ratio ranges from 0.35 to 2; 50% of the map 



O 



is covered by very high ratio, > 1, gas. The very high ratios are predominantly from 
interarm regions and appear to indicate the presence of wide-spread optically thin gas. 
, Star formation tracers are better correlated with the total neutral gas disk than with the 

' molecular gas by itself implying Tisfr ^gas- Using the 100 fim and 21 cm continuum 

from NGC 6946 as star formation tracers, we arrive at a gas consumption timescale of 
^ ' 2.8 Gyr, which is relatively uniform across the disk. The high star formation rate at 

d • the nucleus appears to be due to a large accumulation of molecular gas rather than a 

large increase in the star formation efficiency. The mid-plane gas pressure in the outer 
(R > 10 kpc) HI arms of NGC 6946 is close to the value at the radial limit (10 kpc) 
of our observed CO disk. If the mid-plane gas pressure is a factor for the formation of 
molecular clouds, these outer HI gas arms should contain molecular gas which we do 
not see because they are beyond our detection limit. 

Subject headings: ISM:molecules — stars:formation — galaxies: spiral — galaxies:ISM 
— galaxies:individual(NGC 6946) 
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Introduction 



We have observed N GC 6946 as part of a p r ogram of deep mapping of extended, cold CO 



in nearby spiral galaxies (jCrosthwaite et al 



2001 



20021 ). Large, fully-sampled, and deep images 



of CO in nearby galaxi es are relatively uncom mon. Interferometric images either do not contain 
zero spacing data (e.g., ISakamoto et al.lll999l ) or if they do, the images typically do not go deep 
enough to detect cold, molecular disk gas. The detection of cold, extended molecular gas at large 
galactocentric radii in large spirals is a goal of this mapping project with the NRAO 12 Meter 
Telescope. 



A gas-rich Sc galaxy at a distance of 6 Mpc (jEastman et al.lll996l : ISharina et al. 



2000), NGC 6946 is known for its bright and asymmetric optical spiral arms (jArp 



1997: Karachentsev et al 



19661 ^. The pro- 



nounced a symmetry may be caused by interactions w ith neighboring dwarf galaxies, UGC 11583 



and L149 (Sharina et al 



1997 



Pisano k Wilcotdbood ). N GC 6946 has a 9^ opt ical diameter on the 



sky and an atomic hydrogen (HI) gas disk extending to 25' (iRogstad et al.lll973l ). Like many spirals, 
NGC 6946 has a central ~ 2' diameter "hole" in its HI disk. Radio continuum, far infrared (FIR), 
optical line, and X-ray observations indicate vigorous sta r formatioii in the disk and an interstellar 
medium (ISM) stirred by supemovae and stel l ar winds (|Engargiolal Il99ll : iBoulanger &: Viallefond 
1992 ; Kamphuis fc Sancisi 1993 ; Schlegel 1994 ; Lacev et aP 1997 ) . The high leve l of star formation 
in th e disk of NGC 6946 has been attributed both to its stro ng spiral density wave (ITacconi &: Young 
1990), and to stochastic, self-propagating, star formation ([PeGioia-Eastwood et al.lll984l ). 



Until recently, the CO morphology beyond the inner 3' radius (~5 kpc) was not well known. 
Like many spiral galaxies, NGC 694 6 has a bright nucl e ar CO peak which falls off nearly expo- 
nentiall y with galactocen t ric ra dius (iMorris &: Lol Il978l : lYoung &: Scovilld Il982l ) . A 3' diameter 
map of iTacconi &: YoungI (1198911 re vealed substantial CO emission in the inner disk with several 
emission peaks. ICasoli et al.l (|l990l ) mapped two 2' circular fields in the disk of NGC 6946 in a 
study of arm and interarm regions. The inner 2' region has been mapped with higher resolution 
single dish and aperture synthesis telescopes i n transitions of CO and CO isotopo r ners (iBall et al 



1985 ; IWehachew et al.lll988l : ISofue et al.lll988l : llshizuki et al.lll990l : IWah et al.lll993l : iRegan Vogel 



19951 ). A larger interfero meter-l-single dish map of NGC 6946 has been made by the BIMA SONG 
team (IRegan et al.ll200ll). although w ith lower sensitivity to cold extended emission than the maps 
we present here. I Walsh et al.l (j2002l ) presented 21" resolution fully-sampled single dish maps of 
CO(l-O). We compare our results to theirs. 

The "On The Fly" (OTF) observing mode at the NRAoEl was ideally suited to the imaging of 
extended gas in galaxies. Here, we present 16' by 10' maps of CO(l-O) and C0(2-l) covering the 
optical disk of NGC 6946 made with the 12 Meter Telescope. The maps are deep enough to detect 
cold, extended interarm CO to levels of Ico ~ 1 K km s~^, or ~ 2 x 10'^'' cm~^. We use these 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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primary tracers of the molecular gas phase and combine them with archival and published data in 
order to study the molecular gas disk, the total neutral gas surface density, and their relationship 
to star formation in the disk of NGC 6946. 



2. The Observations and Data Reduction 



2.1. 12 Meter CO Observations 



The observations of CO(l-O) at 115 GHz and C0(2-l) at 230 GHz were made at the NRAO 
12 Meter Telescope at Kitt Peak, on two separate observing runs in December 1996 and March 
1997. An equivalent of ~14 hours on-sou rce observing was a ccumulated for each CO line. Calibra- 
tion was by the chopp er wheel method ( Ulich Sz Haas 19761 ). We convert the recorded T* values 
(jKutner k Ulichlll98ll ) to main bea m temperature, T^T) = T^/Vm^ Vm = 0-88 ± 0.04 at 115 GHz 
and r/^ = 0.56 ± 0.06 at 230 GHz dManguml E996al . Il997l l @. We report Tco = T^b throughout 
this paper. The filter bank spectrometer was configured for 2 MHz channel widths with 256 total 
channels, producing spectral channels of a 5.2 km s~^ width at 115 GHz and 2.6 km s~^ at 230 
GHz. 

The "On The Fly" (OTF) observing mode was used (jManguml Il996bl ') |. Scanning rates 
between 30 and 40 "/sec and scan row spacings of 18" at 115 GHz and 8" at 230 GHz were selected 
to ensure sampling better than Nyquist over the mapped 16' x 10' region. The beamsize (FWHM) 
at 115 GHz is 55" and 27" at 230 GHz. Chopper wheel calibration and sky offs were measured 
every 2 scans. The rms noise level was reduced by averaging the 16 individual OTF maps made at 
115 GHz and 10 individual OTF maps at 230 GHz. 

The OTF data was gridded using the NRAO AIPS package. A linear baseline was removed from 
each spectrum. At this stage of the data reduction the background in the channel maps is uneven 
from one scan row to the next. This spatial variation, due to sky brightness fluctuations between 
successive off-source sky measurements, appears as "stripes" in the continuum level running in the 
scan direction. The maps were made large enough to give us emission free, 2.5' wide, regions at 
the east and west ends of the maps, verified by the examination of individual spectra from these 
regions. These regions were used to remove a linear baseline from each row of the channel maps 
in the scan direction, effectively giving us a uniform background in all the channel maps. The rms 
variation in the channel maps before and after "de-striping" differed by 10 mK at 115 GHz and 30 
mK at 230 GHz; in both cases a 20% reduction in the channel map noise. 



The mean rms noise (la) from fully reduced, line free channels is T^f, = 0.050 K and T^b = 
0.145 K for the CO(l-O) and C0(2-l) respectively. The channels were averaged to a 10.4 km s^^ 



^Available at http://www.tuc.nrao.edu/12meter/obsinfo.html 
''Available at |http:/ /www.tuc.nrao.edu/12meter/obsinfo.html[ 
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channel width; the mean rms noise in the averaged channel maps is 0.035 K for CO(l-O) and 0.085 
K for CO (2-1). The error beam of the 12m Teles cope, which has been measured at 345 GHz, 
is extrapolated to be 8-9' FWHM (jManguml 119971 ). The contributions from the error beam in 
the individual ch annel maps is estirnated to be less than a few percent, which we neglect here 
(see discussion in ICrosthwaite et al.ll2002l ). In the 10.4 km s^^ channel width cubes the rescaled 
uncertainty in T^b due to sky brightness fluctuations is ±7 mK for CO(l-O) and ±11 mK for 
C0(2-l). 

We are interested in bringing out as much of the extended low level CO emission as possible 
in our integrated intensity maps. To do this we want to minimize the noise that would normally be 
added to an Ico map when the channel maps are added together. To test the effects of the level 
of flux cut (level in the individual channels below which the data is deemed to be noise and not 
added to the integrated map) on the final integrated flux, we did the following test. We convolved 
the CO cubes to twice the beam size and used the regions with emission > 3 cr to create mask 
cubes. The mask cubes allow us to define areas of emission using smoothed, hence more sensitive 
images, which is one way of discriminating between signal and noise. We then created Ico maps 
from both the masked and unmasked cubes using flux cuts from 1.0 o" to o" in 0.1 a decrements. 
The total flux in the CO(l-O) integrated intensity map made from the masked cube with a 0.5 a 
flux cut agrees with flux from the unmasked cube with no flux cut, to within 1%. For the C0(2-l) 
map, a more conservative 0.8 o" was used. 

We created maps of the per pixel uncertainty for our Ico maps by combining the per channel 
sky brightness fluctuation and rms noise, counting the number of channels contributing to the 
integrated emission at each pixel, and then adding the uncertainty in r]^ using error analysis 
techniques for the sums and products of independent variables. We created maps of the fractional 
uncertainty in our integrated intensity maps by dividing the aforementioned uncertainty maps by 
the Ico maps themselves. The uncertainty is dominated at the nucleus by the uncertainty in r/^ 
while at the outer edges of the maps the uncertainty is dominated by the rms noise. For the C0(1- 
0) integrated intensity map the combined uncertainty is on the order of 7% at the nucleus, 12% 
along the prominent spiral arm pattern, and 25% for the outer disk emission. For the C0(2-l) map 
the uncertainties are on the order of 11% at the nucleus, 20% along the spiral arm pattern, and 
30% for the outer disk emission. 

To obtain H2 column densities from Icq we applied the "standard conversion factor", Xco 

Xco is estimated to be accurate to a factor 



( Scoville Sanders 



1987 



Young &: Scoville 



1991 



of two in the Milky Way disk (j Solomon et al. 1 119871 . and references therein), although it may 
overpredict m the central arcminute of large spirals (iDahmen et al.l Il998l : iMeier &: Turner 



2OOII ). We assume that these uncertainties hold for NGC 6946 as well. We will adopt the standard 
conversion factor. 



Xco = NhJIco = 2 X 10^" cm"^ (K km s 



-l^-l 



oflStrong et al.l (|l988l ) throughout the remainder of the paper (a recent calibration bv lHunter et al 



19971 . flnds a mean Xco = 1-6 x lO^'^ cm ^ (K km s ^) ^ for the entire Milky Way disk and 
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Xco = 2.7 X 10^*^ cm ^ (K km s ^) ^ for the outer disk, R > R©). Using this conversion fac- 
tor, Iqo can be converted into H2 (no He or heavier elements) surface densities using 

^hMq PC-') = 2.4 X {Ico/ K km s-^) 

including a geometric correction for inclination of the galaxy (40°, Table 1) to correct the surface 
densities to their face-on values. 



2.2. HI, 21 cm Continuum and FIR Maps 



We obtaine d archival VLA D array 21 cm observations of the HI line in NGC 6946, which 
are described in iTacconi &: Younei (|l986l ). The beam size is 49" x 42", with pa = 73°, and each 
channel has a velocity width of 10.3 km s~^. The integrated intensity map was constructed from 
signals greater than 1.7 mJy beam~^ (1.2 a). Absolute flux calibration error is on the order of 
the uncertainty in the flux of 3C 286, < 5%. The single channel rms noise of 1.4 mJy beam~^ 
corresponds to an HI column density of 7.7 x 10^^ 



cm 



Integrated intensities were converted to face-on HI surface densities (no He or heavier elements) 
using, 

J^Hi = 3.3 Mq pc^2 (jy beam-i km s'^)"^ • Ihi- 



From a comparison to the single dish observations of I Gordon et al.l (jl968l ) we estimate 40% of the 
HI flux of NGC 6946 is missing from the VLA maps due to the lack of short baselines (emission 
extended over 15' ). If the missing flux exists as a screen across the disk of NGC 6946 it constitutes 
a missing uniform HI surface density of ~ 1.2 Mq pc~^. 

The total observed HI emissi on in the VLA maps amo unts to Mhi = 5 x 10^ M©, not including 
missing under-sampled emission. iTacconi Sz Youna (|l986l ). using the same VLA HI data, obtained 
3.9 x 10^ Mq, rescaled to D = 6 Mpc. Their mass is probably slight ly lower than ours becau s e they 
used a higher flux cutoff to produce their integrated intensity map. iBoulanger &: Viallefondl (119921 ) 
obtained 6.7 x 10^ (rescaled to 6 Mpc). Boulanger's mass includes single dish observations 
merged with the Westerbork interferometer data to cover the short spacings and corrections for 
Milky Way absorption and emission. The difference between our mass and theirs, 33%, is in line 
with our estimate of the missing flux in the VLA map. 

A Brandt rotation curve was fit to the HI radial velocities, for radii < 10'. Our resulting 
kinematic paranieters , Table 1, are in good agreement with the rotation parameters found by 
Carignan et al.l (|l990l ). A total dynamical mass estimate based on the fit (see Table 1) is: Mdyn = 
1.9(±0.6) X 10^^ Mq. 



^The adopted inclination was derived from a Brandt model fit to the HI data. 
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We produced a 21 cm continuum map by averaging 4 "line free" channels in the NGC 6946 HI 
channel cube. The la noise level is 1.2 mJy beam~^. 

High resolution (HiRes) 60 and 100 //m IRAS maps were obtained from IPAC 0. The 100 //m 
HiRes map was produced by 200 iterations of a maximum entropy method yielding a 69" x 65", 
pa = 20° beam and a la noise level of 2.0 MJy sr'^ The 60 /im HiRes map was produced by 50 
iterations of a maximum entropy method yielding a 45" x 41", pa = 21° beam and a la noise level 
of 0.6 MJy sr~i. 



3. Images of CO in NGC 6946: The Molecular Gas 



A relatively heavy molecule, CO is easily excited at the temperatures of molecular clouds, 
with its first rotational level at an energy of E/k = 5.5 K. Since it has a small dipole moment 
and is generally optically thick, the excitation of the CO rotational levels is driven by collisions 
at H2 densities above Ucrit ~ 300 cm~^ for the J=l-0 transition. The J = 2-1 transition, with an 
upper level energy of E/k = 16.6 K, has a higher critical density, although its higher optical depth 
somewhat compensates for that. For gas kinetic temperatures of ~7-10 K, we expect C0(2-l) to 
be thermalized at critical densities above ricrit ~ 10^ cm~^, slightly higher than CO(l-O) values. 
One might therefore expect C0(2-l) emission to be more tightly confined to the spiral arms than 
CO(l-O), since the arr ns probably have den ser gas, and are also closer to HII regions that might 
warm the clouds (e.g., IScoville et al. 1119871 ). However, as we shall see, the C0(2-l) and CO(l-O) 
properties in NGC 6946 are not precisely what one might expect based on Galactic clouds. 

We do not present CO kinematic data (radial velocity maps, veloci ty dispersioii maps , rotation 
curves). Our kinematic results do not differ significantly from those of I Walsh et al.l (|2002l ). and we 
refer the reader to their presentation of the kinematics of NGC 6946 derived from IRAM single 
dish data. 



3.1. CO(l-O) Maps 

The 24 CO(l-O) channels containing line emission are displayed in Figure [H The channels 
have been smoothed to a 10.4 km s~^ channel width. The channel maps show the butterfly pattern 
characteristic of emission from an inclined rotating disk. High velocity nuclear gas shows up as 
emission at the location of the nucleus spread across nearly the entire bandwidth of the line emission. 

The disk of CO emission in NGC 6946 extends 8' (14 kpc) north-south by 10' (18 kpc) east- west 
as shown in the integrated intensity, /io= / Tiq dv (Figure [2^), and the peak main beam brightness 



^Descriptions of the IRAS HiRes data reduction and HiRes data products are available at 
|http://www.ipac.caltech.edu/ipac/iras/toc.html[ 
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temperature, Tiq, (Figure [2)3) maps. I Walsh et alj (|2002l ) also found CO emission extended out to 
10' in 22" resolution single dish IRAM observations. There are no significant difi'erences between our 
and Walsh's map that cannot be attributed to the difi'erence in resolution. The CO(l-O) emission 
covers the entire optical disk of a Digital Sky Survey (DSS) image of NGC 6946 (Figure [J]:) . 

The total CO(l-O) fiux density measured over a velocity width of 250 km s^^ is 1.2(ib0.3) x 
10"^ Jy km s~^, in agreement with the value published in the FCRAO survey of 1.2(ib0.4) x 10^ 
Jy km s~^ (jYoung et al.lll995l ). but hi gher than the 0.7 9(=fc0.01) x 10^ Jy km s~^ detected in a 



22" resolution single dish IRAM map (jWalsh et al.ll2002l ). Converting our total observed Iiq into 



a total molecular mass, we obtain M//2 = 3.3 x 10^ Mq for NGC 6946; the mass will be higher 
if there is outer disk C O that we do not detect. This is consistent with the mass obtained by 
Tacconi fc Younel (|l989l ). 3.0 x 10^ Mq (rescaled to 6 Mpc and our Xco)- 



The CO integrated intensity (/lo) and temperature maps (Tiq) both show the same features: a 
bright nucleus, a broad asymmetrical spiral arm pattern and the overall asymmetrical distribution 
of the disk emission. The nucleus has a CO emission peak at a = 20^^ 34^52. 6* 6 = 60°9'l6" 



(J2000), with Iio"""^ = 72 K km s'^, TV^a""'' = 1-4 x 10^2 cm'^, or in terms of surface density. 



170 Mq pc (face-on value), averaged over our 1.6 kpc d i amete r beam. The total 



molecular mass within the central beam is 4.7 x 10^ Mq. Ilshizuki et al.l (1990) obtained 4.5 x 10*^ 
Mq (adjusted to 6 Mpc and our conversion factor) in a 65" diameter region from interferometric 
observations that recovered 70% of the total flux. Large scale CO spiral structure is traced in a 
broad "S" pattern, outlined by the 7.5 K km s~^ contour in the Iiq map and more clearly seen in 
the Tio map, outlined by the 0.23 K contours. Here, molecular gas surface density is > 20 Mq pc"^. 
There is also extended, weaker CO emission from the disk outside of the 7.5 K km s~^ contour in 
the IiQ map. Large regions of molecular gas exist to the northwest and southeast of the nucleus, 
as well as patchy emission out to the edges of the maps. Not all of this gas is associated with the 
weaker optical arm pattern. The face-on gas surface density for this extended region has a mean 
Ico of IiQ ~ 3 K km s~^ corresponding to Nhj = 5 x 10^'' cm"^ or E/^j ~ 10 Mq pc"^. This 
extended portion of the gas disk contributes 1/4 of the total CO emission we observe. The total 
mass of this extended gas disk is Mh^ ~ 10^ Mq. 



3.2. CO(2-l) Maps 

The 24 channels containing C0(2-l) emission are displayed in Figure El The rms of these 10.4 
km s~^ channels is three times as high as for CO(l-O), so the maps are less sensitive to faint emission. 
The total C0(2-l) flux density measured over a velocity width of 250 km s'^ is 3.5(±1.0) x lO'' 
Jy km s~^. The butterfly patterns for the C0(2-l) and CO(l-O) are generally the same, with the 
differing resolutions, 27" and 55" respectively, accounting for most of the differences. 

As with CO (1-0) the bright nucleus and asymmetric distribution of arm and disk emission are 
also apparent in the integrated intensity, I21, and peak main beam brightness temperature, T21, 
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maps (Figures andllja). The C0(2-l) emission in NGC 6946 is almost as extensive as that of 
the CO(l-O): the north-south and east-west dimensions are only 1' shorter. The observed C0(2-l) 
emission covers the brighter portions of the DSS optical image of NGC 6946, Figure Hh. 

The nuclear peak oftheCO(2-l) integrated intensity map is at a = 20''34'"53.1^ 6 = 60°9'l5" 
(J2000) with l2i"'"'^ = 125 K km s~^. Using a Galactic disk Xco^ this corresponds to a peak col- 
umn density of Nh2^^^ = 2.2 x 10^^ cm~^ or 450 Mq pc~^ averaged over the 790 pc beam. If other 
nearby galaxies are any indication, this may be an overestimate of the true molecular gas mass in 
the nucleus, by factors of ~ 2 — 3(§3.1.) The largest T21 value, 0.85 K, occurs at the same location 
as the I21 peak and both peaks are consistent with the location of the /lo peak. 

At the resolution presented of maps, inside a 2' (3.5 kpc) radius both the C0(2-l) and CO(l-O) 
emission trace a broad 2 arm pattern. At a 2.5' radius (4.4 kpc) the C0(2-l) still traces a broad 
2 arm pattern while the CO (1-0) emission from a larger beam has clearly bifurcated into a 4 arm 
spiral (Figure [5^). At a 3.5' radius (6.1 kpc) 5 arms are discernible in azimuthal emission plots of 
both CO emission lines (Figure [5)3). 

Along the optical arm pattern the C0(2-l) emission is uneven, (Figure Hh) with typical varia- 
tions of 15 and 5 K km s~^ separated by 40" or less. This same c lumpy patter n of C O arm emission 



can be seen in the 22" CO(l-O) and CO(3-2) maps presented by I Walsh et al.l (120021 ). The emission 



on each arm at the same galacto-centric radius is also uneven as can be seen in Figures [5^ and b. 

At 27" resolution it becomes possible to distinguish between arm and interarm emission. The 
I21 contrast between the brighter, on arm, clumps and obvious off-arm regions is ~ 3. Using the DSS 
optical image to mask the outer (R> 2.5') arm/interarm regions, we obtain a mean arm/interarm 
contrast of 2.4. This i s higher than both the 1.2 and 1.8 values found from CO(l-O) and CO(3-2) 



TRAM observations of I Walsh et al.l (120021 ) who interpret the increase in contrast at the higher CO 



transitions verses CO(l-O) as evidence that the on-arm molecular gas is warmer. 

The presence of widespread interarm gas becomes plainly apparent in a comparison to an opti- 
cal image, Figure U^. The fact that interarm C0(2-l) is present and at emission levels comparable 
to CO (1-0) suggests that the interarm gas is also not very cold and that i t must be warmer than 



~ 4 — 7 K. This widespread distribution of interarm gas is not seen in the lSakamoto et al.l (119971 ) 
C0(2-l) survey of the Milky Way. 



3.3. Excitation of CO in NGC 6946 

A CO ratio map, ri2 = l2i/IiQ^ is presented in Figure [6^. The C0(2-l) channel cube was 
convolved to the beam size of the CO(l-O) data and the > 2>(j emission from the beam matched 
cubes was used to produce the ri2 map. The spiral pattern is only weakly reflected in the ri2 map. 
Outside of the inner few arcminute radius the pattern breaks up into a patchwork of high and low 
ri2 regions. 
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The ratio uncertainty for the regions of strong emission is dominated by the error in the 
conversion to main beam temperature, while the uncertainty for the weaker emission regions is 
dominated by the rms noise in the contributing channels of emission. Adding in quadrature the 
fractional uncertainties from the Iqo maps we obtain an error estimate on the order of 15% at the 
nucleus, 25% along the spiral arm pattern, and 40% for the outer disk ratios. We were concerned 
that the ratio might be biased towards higher values due to the uncertainty in the denominator. 
We performed a standard statistical analysis to correct a ratio for bias this bias. The bias corrected 
ri2 map was not significantly different from the map presented in Figure Eti and well within our 
reported error. Despite the large unc ertainty estimates, our ratios are consistent with those derived 



from other independent observations (Weliachew et al 



Wall et al 



1993! ^. 



1988 



Ishizuki et al 



199C 



Casoli et al. 



199C 



In principle, ri2, the ratio of C0(2-l) to CO(l-O) integrated intensities, gives information on 
the excitation of the CO emission and depending on circumstances, the kinetic temperature, T^. 
The observed main beam temperature at frequency, is related to a Planck function, B^^ with 
excitation temperature, T^x, by the areal beam filling factor, f^, a conversion to the Rayleigh- Jeans 
temperature scale, and radiative transfer through a optical depth, Ty-. 



f mb 



including a correction for the cosmic microwave background contribution to the beam, By{Tcmb), 
at temperature 2.73 K. Standard assumptions are that f^ is the same for both emission lines, that 
both lines are optically thick (r,y ^ 1) and that the emission ca n be characterized by a single 
Tex (IDickman et al.l Il986l : iMaloney &: Blackl Il988l : ISakamotd Il996l ) . With these assumptions we 
can naively derive T^x values for the emission: ri2 = 0.9 corresponds to T^x = 20 K, while ri2 = 
0.5 corresponds to T^x = 3-5 K. Higher ratio values are found in the nuclear regions of starburst 
galaxies where elevated star formation heats the molecular gas; lower values are typical of cold or 
subthermally excited (Tex< Tfe) disk molecular clouds. Values of ri2 > 1 indicate optically thin gas 
in this picture. 



Within 1' of the nucleus, ri2 ^ 0.8. This is consistent with other studies (jWeliachew et al 



1988 



Ishizuki et al 



199C 



Wall et al.lll993l ) that indicate the presence of relatively dense molecular 
gas a s temperatures of abo ut Tex ~ 10 — 15 K. This is comparable to the Milky Way mean ri2 of 



0.74 (|Qka et al.lll996l . Il998l l 



Along the optical arms the mean ri2 is 0.83, ranging from 0.44 to 1.4. To within our uncer- 
tainties, this is the same as the 0.73 mean found for M ilky Way molecular arms where a similar 
range of ri2 values is found ([Sakamoto et al.lll995l . 119971 ). Prom a naive interpretation of the ratios, 
a mean Tex of 12 K is implied for the arm CO in NGC 6946. 



Casoli et al.l (|l99d ) concluded that the interarm ri2 does not differ significantly from that of the 



arms based on observations of two circular regions in the arms, one on the east side of NGC 6946, 
the other on the west side. The two circular regions they observed have a mean ri2 of 0.8 in our 
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ratio map, also not significantly different from that of the arm regions. However when larger regions 
of the CO disk are sampled, a very different picture emerges. There are regions (5% of the ratio 
map) where ri2 < 0.6 indicating cold {T^x < 5 K) or sub-thermally excited CO. There are also 
much larger regions where ri2 > 1. 

Very high ratio, ri2 > 1, gas covers 50% of our ratio map; and these are predominantly 
interarm regions. Figure [6)3. We have eliminated the possibility that these very high ratios are 
due to contributions from the telescope error beam at 230 GHz: the error beam at the 12 Meter 
Telescope is a measured quantity and modeling of the error beam contribution in the individual 
channel maps shows that it contributes at most a few percent of the emission in the I21 map. 

The high ri2 ratios appear to indicate the presence of widespread, warm, optically thin CO- 
emitting gas in the interarm regions. While optically thin CO emission is not typical of Galactic 
molecular clouds, our large beam and unusual perspective from outside the disk of NGC 6946 may 
allow us to trace a component that is difficult to detect in the Galaxy. We note that the photon- 
dominated regions (P DR) around cool stars can produce bright CO emission with high temperatures 
(|Spaans et al. Ill994l ) since their sof ter radiation fie l ds all ow the CO to be photoelectrically-warmed 
at low A^. Also, as pointed out bv lWiklind et al.l (|l99d ) it does not take much optically thin gas 
to swing the emissivity; the C0(2-l) emission that we see could easily come from a small amount 
of optically thin, but "CO-l oud" gas that is not the dominant cloud population. Other less likely 
possibilities are discussed in ICrosthwaite et al.l (j2002l ) . 



We have also seen a have a high overall ri2 ~ 1 in M83 (ICrosthwaite et al.1 120021'): confirmin g 



earlier reports of unusually high ri2 in the se two galaxies (ICastets et al 



Larg e scale ri2 maps of the Milky Way (jSakamoto et al 



199 



I 



1990 



Oka et al 



Wiklind et al.lll99dl 



19961 : 



Sakamoto et al, 



19971 : lOka et al.lll998l ) show high ratios in immediate proximity to pockets of star formation. The 
widespread regions of high ri2 seen in M83 and NGC 6946 are not seen in surveys of CO in our 
Galaxy. 



4. CO, HI and the Neutral Gas Surface Density 



Although CO emission covers only 15% of the HI disk surface area, NGC 6946 is nevertheless 
a molecular gas-rich gal axy; 1/3 of the hyd rogen gas disk is in molecular form. This is consistent 
with the 36% value from lWalsh et al.l (j2002l ) once the flux missed in the interferometer HI maps is 
taken into account. 

In Figure [7^, we present our CO(l-O) emission in contours superimposed on a grey scale image 
of the HI emission. This image shows that CO emission fills the inner 6' of the HI disk. The HI 
disk has a central depression, which rises to peaks in a clumpy ring of emission. The HI ring is 
coincident wit h the outer boundary of CO(l-O). While spiral structure is readily apparent in the 
outer HI disk (jTacconi &: Young 119861 ) . it is insignificant in the inner disk, R< 5' (9 kpc), where 
the HI column densities decline and H2 takes over. At 1' resolution there are no obvious systematic 
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displacements of the HI relative to the CO(l-O) spiral structure other than the strongest HI tends 
to be near the outermost ends of the CO(l-O) arms. In fact, if anything, HI and CO peaks tend to 
coincide, which is more apparent in the higher resolution C0(2-l) maps. 

The CO (2-1) image in contours is shown in Figure [Zb, again with a grey scale image of HI. 
Outside the inner 2' (3.5 kpc) region of the HI depression, the peaks of C0(2-l) emission tend to 
coincide with HI peaks along the ring. This could be evidence for photod issoci ated H2 creating a 
higher HI column density at these locations. iTilanus &: AllenI (|l989l . Il993l ) and lWilson &: Scoville 
(|l99ll ) suggest the same explanation for HI in the spiral arms of M51, M83 and M33. 



From these images, we have constructed plots of the azimuthally-averaged gas surface densities, 
Sj^2 5 ^Hi and ^gas, as a fu nction of galactocentric radius in Figure [8] (corrected for inclination). 
This is similar to the plot of Young &: Scovill3 (19821), w hich was based on a radial cross map, and 



those of lTacconi fc Youn j (|l98d ) and 



Walsh et al 



(|2002l ) . H2 dominates the neutral gas distribution 



in within the central 3' (5 kpc) of NGC 6946, rapidly fading into a predominantly HI disk beyond. 
A least squares fit of an exponential to Ic q ove r the entire r ange of radii has a scale length of 
1.0' ±0.1 (1.8 kpc). iTacconi Youn j (|l986l l and lWalsh etaP ^200^ ) obtained scale lengths of 1.5 ' 
and 1.4' for the CO disk from FCRAO and IRAM single disk observations; iRegan et al.l (|200ll ) 
obtained 1.3' from the BIMA interferometer survey data combined with NRAO 12m observations. 
All these latter scale lengths exclude the nuclear contribution to the fit and are sensitive to the 
ra nge of radii over which th ey co mputed. We obtain 1.2' and 1.3' fitting to the range of radii used 
bv iTacconi fc Younsl jl98g) and Iwalsh et al.l (|2002l ). 40" - 280" and 60" - 260" respectively We 
conclude that the CO falls off with an exponential scale length of 1.2' it 0.2' on arcminute scales, 
excluding the inner 1' nuclear region and about 1' if the nucleus is included. 

A total neutral gas surface density map, Sgas= 1-36 {T,hi+ is presented in Figure [9^. 

This map has been corrected for inclination but not for missing 'Shi due to under-sampling in the 
VLA image (~1 Mq pc~^), A false color image is presented in Figure [9)3, with HI in red and CO 
in green, so that the molecular contribution to the total gas disk can be distinguished. Values for 
T,gas range from 240 Mq pc~^ at the nucleus, to 20-45 pc~^ along the spiral arms within the 
optical part of the galaxy, to 5-10 Mq pc~^ for the outer arms. At first glance, H2 and HI blend 
into a single global gas structure extending from the nucleus to the outer gas disk. A closer look 
reveals some subtleties. In the southern part of NGC 6946, the neutral density gradient is smooth 
and gradual. In the north, the gradient is much steeper, with the gas density falling to a level of 
20 Mq pc~^ in 90" (2.6 kpc) to the north, into an interarm void, as compared to a 180" fall to 
the same mass density in the south (this can also be seen in Figures [Tj). This steeper falloff on the 
northern side of the central gas peak may be related to the more well-define d spiral arm to the 
north, and perhaps also to the prominent northeastern optical arm (lArpl ll96d). 



The global properties of three late galaxies for which we have a similar data set, along with the 
Milky Way for reference are listed in Table 2. IC 342, M83 and NGC 6946 all have gas disks that 
extend well beyond their optical diameters and all three show some evidence for a bar. The Yjh^'^^^^ 
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values for IC 342, NGC 6946 and M83 are all significantly higher than that of the Milky Way; a 
factor of four higher for IC 342 and NGC 6946, a factor of seven higher for M83. The NGC 6946 
Y,^^disk ^^^^g similar to that of IC 342, while its E^a'^'"' intermediate between IC 342 and M83. 
All four galaxies have similar relative proportions of H2 an d HI, global M^^/Mfff values between 
0.3 and 0.5 are also seen in the nearby galaxy, Maffei 2 (jMason &: Wilson 1 12004| ). M83 stands 
out as a particularly gas-rich object relative to its dynamical mass, with a factor of two higher 
Mgas/Mdyn than IC 342 and NGC 6946, and four times that of the Milky Way The II2 and optical 
disk diameters of the three external galaxies are the same, Dij2/D25 ~ 1- The dynamical masses of 
these three galaxies are an order of magnitude smaller than that of the Milky Way, although their 
total gas mass is similar. It is difficult to place too much emphasis on subtle variations that rely 
on distance estimates, since the uncertain distances to these nearby galaxies can easily affect their 
masses, both gaseous (oc R'^) and dynamical (cx R). 



5. Star Formation and the Neutral Gas Disk of NGC 6946 

5.1. Star Formation Tracers, CO, and the Neutral Gas Disk 

Because stars form from molecular gas, we expect to see correlations between CO and tracers 
of star formation. These tracers include: Ha emission from ionized gas associated with massive 
young stars, FIR emission from dust heated by young stars, and radio continuum emission from 
relativistic electrons created by supernovae. To calculate linear correlation coefficients all the maps 
were convolved to the largest beam size in our collection of maps, a 70" circular beam size, then 
dipped to use > 3a emission, and sampled on a 35" grid. 



Ha emission (jFerguson et al.lll998l ) and I21 are compared in Figure [TOk. The nuclear bar traced 
by Ha within the central 2' is aligned with the elongation of the central I21 contours. The I21 arms 
are traced at higher resolution by the Ha emission. In the outer I21 disk, the brighter patches of 
C0(2-l) emission are aligned with the brightest patches of Ha which is most clearly seen along the 
northeastern spiral arms. This supports the conjecture that the ring of HI emission seen at these 
radii is the result of photodissociated H2. It is apparent from this map that ratio of Ha emission 
to I21 is not the same in the nucleus and outer arms; due to substantial extinction particularly in 
the nucleus. Ha is probably not the optimal tracer of star formation for this gas-rich galaxy. When 
Ha is compared to I21 or our primary molecular mass tracer, Iiq, we obtain a linear correlation 
coefficient of ~ 0.65. The Ha correlations do not improve when a 1' diameter region centered on 
the nucleus is excluded from the correlation. 

Figure llOb shows that the 60 //m FIR emission generally reproduces the morphology seen 
in I21 (grey scale), despite the noticeable artifacts (the "boxy" structure of the centr a.1 contours) 



produced by the HiRes maximum entropy algorithm. An ISOPHOT 60 /um map by iTuffs et al 



(jl996l ) confirms the overall morphological structure seen in the HiRes map. The 60 /xm nuclear 
peak is at a = 20''34™52.4'' 5 = 60°9'l5" (J2000) which is, within the uncertainties, coincident 
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with the CO peak. If C0(2-l) along with 60 fim FIR trace a warmer component both heated 
by star formation; we expect these emission pairs to have similar feature in the disk. The 60 /im 
emission traces the overall C0(2-l) arm structure better than the 100 nm emission shown in Figure 
[TOb . The 4 prominent peaks in the outer 60 /im emission disk are aligned with peaks in the outer 
/21 arm structure. 

The 100 fim FIR and 21 cm continuum maps have resolution comparable to our CO map 
(Figures [TOb and d). The 100 //m emission roughly mimics the CO(l-O) emission disk, but the 
CO arm structure east of the nucleus is not well traced by the 100 fim. The 100 /um nuclear peak 
is at a = 20'^34'"53.0^ 6 = 60°9'23" (J2000); northeast of the CO emission peaks. The radio 
continuum does a better job reproducing the CO morphology, with the exception of the 21 cm 
peak 4' northwest of the nucleus. This is consistent with the suggestion that cosmic rays may 
be more important for the heating of molecular gas on large scales such a s observed here tha n 



the radiation from in situ massive star formation (Adler. Allen. &: Lo 



1991 



Suchkov et al.lll993l ). 



The continuum peak to the north and w est of the nucleus at a = 20'^34™24.9'' 6 = 60° 10 38 



is probably a background radio galaxy (Ivan der Kruit et al.l Il977l . based on its spectral index, 
a = -0.84 and flux.) The 21 cm continuum nuclear peak at a = 20^* 34^52. 4* 6 = 60°9'l5" is 
consistent with the CO(l-O) peak. Much of the outer CO arm morphology and prominent emission 
peaks are traced by non-thermal radio continuum. 



Walsh et al.l (120021 ) also found statistically significant correlations between CO(l-O) emission 



and the total radio-continuum at 6 cm, and between CO(l-O) and 20 cm no n-thermal cont inuum 
in NGC 6946. When fluxes are extrapolated to the expected 20 cm values (ICondonlll992l . using 
Si, oc z^"'''^), we obtain a mean /10/T20 of 1.4 ± 0.5 km s^^ from our 70" beam convolved maps for 
the di sk of NGC 6946. This is higher that the /10/T20 =1-0 value determined bv lAdler. Allen. &: Lo 
(jl99ll ) from sparely sampled Jio data, but in line with the 1.3 mean value for all the galaxies in 
their sample. 

One might naively think that the FIR has a more direct causal link to newly formed stars, 
while radio continuum emission, which arises from relativistic electrons produced in supernova 
remnants, would be more spatially removed. In reality the far-infrared cirrus and non-thermal 
disk emission on arcminute scales represent a less temporally localized reflection of star formation; 
the close correspondence of radio continuum and CO supports a model of cosmic ray confinement 
on 1 kpc scale lengths, and highlights the importance of gas surface density in coi itaining the 



magnetic field required to extract synchrotron emission from relativistic electrons (jBicav et al 



1989 



Bicav fc Heloulll99d : iHelou fc Bicavlll993l ). 



Table 3 lists least squares fits to the radial distribution of ratios of Iiq or I21 to our star 
formation tracers, Gastracer /SF tracer ■ With the exception of Ha, all the star formation tracers 
are highly correlated with CO emission. From the 70" beam maps we obtain linear correlation 
coefficients on the order of ~ 0.95 or better in comparisons between either of our CO maps and 
any of the star formation tracers except for Hq. The radial gradient is reduced by a factor of two 
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when IiQ is used instead of I21 indicating Iiq remains the best tracer of molecular gas associated 
with star formation on kpc size scales. 

While CO is highly correlated with star formation tracers, the correlation of total gas, HI+CO, 
is even better (see Table 3). The radial gradient is reduced by a factor of ten when T,gas is used. If 
we assume that the 21 cm and FIR emission trace star formation, this result implies that the star 
formation rate in NGC 6946 is more closely tied to J^gg^ than S fj-, ■ R estating this result: a simple 
Schmidt law applies, SFE oc T.gas"' (jSchmidtl Il959l : iKennicuttI Il998bl ) . at least on kpc size-scales. 
Moreover, it appears that in the outer parts of NGC 6946, star formation is likely to take place in 
regions of predominantly HI gas. 



We obtained similar resuhs for IC 342 and M83 (ICrosthwaite et al.ll200ll . l2002l ). While the 



correlation coefficient for Y^gas /S21 cm ill M83 is also high, 0.95, for IC 342 the correlation coefficient 
is only 0.73. The lower map correlation for IC 342 is related to the peculiar 1 armed appearance of 
the 21 cm continuum map that contrasts starkly with its grand design appearance in other tracers. 
From least squares fits to the Ygas/^2icm ratio we obtain slopes of —0.012 ±0.004 (dex) for IC 342 
and -0.13 ± 0.02 (dex) for M83. For M83, the slope is not as flat as in IC 342 and NGC 6946. 
This gradient in M83 may indicate an age difference with large r T,gas/S2icm values resulting from 
younger star forming regions in its disk (ISuchkov et al.lll993l ). The Sgas/S2icm comparison for 
NGC 6946 suffers none of these peculiarities. 

A complication in the interpretation of molecular gas mass could arise because of a radial 
gradient in metallicity an d because the Galactic CO conversion factor over-predicts the H2 mass in 
the nucleus o f NGC 6946 (IMeier k Turneilbood ). Based on the metalhcity gradient for NGC 6946 
me asured by iBellv &: Rovl (|l992l ) and the metallicity adjustment to the CO conversion factor used 
by IWilsonI (| 199 il l for giant molecular clouds in M33, Xco could be a factor of two higher at the 
edge of the NGC 6946 CO disk. This would increase the molecular surface density derived from 
CO at the disk edge by a factor of two, and potentially reduce the radial gradient of molecular 
mass implied by Iiq with radius. However, this variation in Xcn does not appear to be universal in 
low metallicity systems (IMeier. Turner. &: Beckll2002l : iLeroy et al.ll2005l ). so we do not adopt the 
correction here. 



5.2. The Star Formation Efficiency of the Disk of NGC 6946 



Studies of the links between star formation and gas, or star formation efficiency, have recently 
focused on corre l ations as a function of position in the galaxy, as opposed to global measures. 



Rownd &: Yound (1199911 examined the relation between T,h2 ^-rid the star formation rate, SFR, as 



traced by Ha; iMurgia et al.l ([2002) compare a radio continuum based S FR a. nd T1H2 '1 IWong &: Blitz 



([2002) compare a SFR from Ha emission to T,gas- lYoung &: Scovilld (|l982l ) pointed out that the 
falloff of T,H2 in NGC 6946 with radius was very similar to the falloff in B luminosity, implying a 
constant star formation rate per unit H2. The most widely used formulation is the Schmidt law. 
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Tigas (|Schmidtlll959l ): and the most widely accepted calibration was done bv lKennicuttI 
(|l998al l who finds A = 2.5 ± 0.7 x lO""^ and iV = 1.4 ± 0.15 for T^gas in M© pc'^ and Y^sfr in 
Mq yr~^ kpc~^. How well does the Schmidt law describe star formation in NGC 6946? 

We can estimate star formation rates from our FIR and 21 cm continuum data. Unlike Ha 
emission, these tracers are relatively unaffected by absorption. Although cirrus heated by an older 
population o f stars may contarninate FIR in some galaxies, this is probably not a big effect in Scd 
type spirals (jSauvage &: Thuanlll992l ). The 21 cm continuum appears to be a good star formation 
tracer on global scales, although it must break down at smaller scales. The maps for this analysis 
were convolved to the the largest circular beam in the set, 70". 

We built a FIR map from the 60 /um and 100 ^m IRAS maps using: 

Sfir (W m-2 sr-1 ) = 12.6 x lO'^^ (2.58 S^o^m (Jy sr"!) + Sioo^m (Jy sr-^)) 



(jCatalogued Galaxies and Quasars in the IRAS Surveylll985l ). We adopt the FIR star formation 
rate: 

SFR^,, (Mo yr-i) = 4.5 x 10-^7 Ly,,(W) 



(lKennicuttlll998bl '). where L = AttD^S and D is the distance to NGC 6946. We derive the T,sFRfir 
map from the Sjir map using: 



^sfRjJMq yv-' pc-') = 5.38 x 10"^ Sf„ (W m"^ sr"^ ) cos i 
where cos{i = 40°) gives us face-on surface values. For the non-thermal continuum at 21 cm we use 



SFR^te (M© yr-i) = 1.9 x 10^ 



-22 ^0, 



(GHz) LntciW Hz-i) 



to get the star formation rate (|Condonlll992l ) and use: 

^SFR^tS^e yr"^ pc-2) = 2.30 x 10"^ SiA2Ghz (Jy beam"!) cos i 
to build the T,sFR„tc map. The T,sfr maps were azimuthally averaged in 30" radius bins. 

The results are show n in Figure [TTk . The SFRs derived from the FIR and radio continuum 
differ by a factor of two. iBelll (j2003l ) finds that the non-thermal continuum from low luminosity 
galaxies is substantially suppressed, and so calibrations of the SFR from radio continuum from 
surveys that do not take into account the overall luminosities of the constituent galaxies tend to 
underestimate the SFR. Bell provides an alternate calibration: 



SFRntc (Mq yr-1) 



5.5 X 10^22 
0.1 + 0.9(L/L,)0-3 



Lntc (W Hz-i) 



where the total galactic flux satisfies the condition, L < Lc, Lc = 6.4 x 10^^ W Hz ^. In the case 
of NGC 6946, L = 5.2 x lO^^ W Hz^i. This cahbration of SFR„te is also shown in Figure EJi 
and agrees with the FIR derived rate. Although the agreement with the SFRjjj. result is clearly 
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promising, this alternate calibration was established using total galaxy fluxes and has not been 
studied with respect to the varyi ng non-the r mal co ntinuum fluxes within galaxies. Figure [TTb also 
shows the SFR predicted by the iKennicuttI (|l998al ) calibration of the Schmidt law using our Sgas 
map; in NGC 6946, it overestimates the SFR relative the SFRjjj. and SFR„fc values. 

Another way to look at the Schmidt law is in terms of the star formation efficiency, SFE. In 
Figure [TTb we express the SFE for each of our SFRs as tsfe, the time it would take to consume 
all the gas at the given star formation rate, tsfe = ^SFRf^gas- Using the SFRj^^ and Bell's 
calibration for SFR^tc, the mean tsfe for NGC 6946 is 2.8 ± 0.8 Gyr and S.s fr oc S^fj,.s. This 



isconsistent with the mean tsfe of 2.6 Gyr for undisturbed galaxies found bv iRownd &: Young 
(jl999l ). The = 1.4 Schmidt law is a lso shown on Figure [TTb : the Schmidt law over-predicts 



the SFE and its dependenc e on S^n.^ . IWong Sz Blita (120021 ) found N values as low as 1.1 for 



i I — 1 ■ I ' rii = <^ — — — •/ 

the galaxies in their sample. lElmegreenI ( 20021 . see references within) has noted that the Schmidt 



law is inconsistent with observations that the SFE is roughly constant across range of galactic 
environments. This certainly is the case for NGC 6946. 

On kpc scales, the SFE for the nucleus of NGC 6946 does not differ from that of the disk; 
in other words, the strong central FIR emission from the nucleus NGC 6946 results from having 
more molecular material from which to form stars rather than any increase in the star formation 
efficiency. This can be seen in Figure [12] where we show the radial distribution of the ratio of S^as to 
21 cm continuum; the ratio does not decline at the nucleus. The lowest values (highest efficiencies) 
are found between 2' and 4' radii (3.5 to 7 kpc). The highest values (lowest efficiencies) are at the 
largest radii. In two other late type galaxies for which we have a similar data set, IC 342 and M83, 
there is a noticeable drop in log{'Egas/S2icm) from the peak values in the disk to the nucleus; an 



indic ation of a elevated star f ormation efficiency relat ive to the disk rates (jCrosthwaite et al. 



2001 



2003). Using the definition of lRownd Yound (jl999l ). NGC 6946 would not be a starburst galaxy 
because it does not have a gas cycling time shorter than 10^ years. Other authors have com e to a 
similar conclusion about the star formation efficiency in NGC 6946: [Meier Turned (|2004l ) based 
on a determination of the nuclear molecular mas s froni high resolution observations of CO isotopes 
and multiple line transitions, and [Madden et al.l (jl993[ ) by evaluating the SFR using CII maps. 



5.3. Is There an Edge to the CO Disk in NGC 6946? 

Does the outer edge of our CO(l-O) map represent a real threshold for the formation of molec- 
ular clouds? Clearly in NGC 6946, there is not a steep or marked decline in the gas as we detected 
i n M83; that galaxy ha, s a real gas edge, in both CO and HI, possibly due to tidal interaction 
(jCrosthwaite et al.[ [20021 ). In NGC 6946 the exponential disk continues as far as we can detect it 
(Figure [8|) . But there may be other subtle changes in the gas with radius and we examine these 
below. 



We use our Galaxy as guide on what to expect for the radial falloff in gas and star formation 
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in a spiral galaxy. The number of HII re gions per kpc^, ^h-, , and at i?Q and 12 kpc (rescaled 
to Rq = 8 kpc) for the Milky Way from IWouterloot et alj (|l988l ) are presented in Table 4. In the 
Galaxy, between 8 and 12 kpc Nhii /^pc^ has fallen off by a factor of 15, Sj^j by a factor of 20, 
while TiHi sta ys relatively fiat|^ Star formation per unit H2 remains relatively constant over this 
range of r adii. IWouterloot et al.1 (jl988l ) found evidence for molecular outflows at distances as large 
as 16 kpc. iLequeux et al.l (119931 ) detected Galactic molecular gas in a search for CO absorption lines 
against continuum sources; their findings suggest cold molecular gas may be far more abundant at 
large radii than previously suspected; t hey suggest the molecular gas mass may be 4 times the HI 
mass at 12 kpc. iMead Kutnerl (|l988l ) detect a sparse distribution of Milky Way molecular clouds 
at 13 kpc with a mean T* ~ 3 K and a mean radius of 20 pc. 

The edge of the observed CO disk in NGC 6946 lies midway between 8 and 12 kpc for an 
assumed distance of 6 Mpc. Comparing NGC 6946 to the Milky Way (Table 4): at 8 kpc (4.6' ) 
Nhii/^VC^ is a factor of four less in NGC 6946 than the Galactic value at Rq, and could be due to 
extinction in Ha. is same for both galaxies (within the errors), and S/// is at least twice 
the Galactic value. A t 12 kpc, Nhtt l^vc^ is the same as Galactic (extrapolating 20" in radius past 
the radial limit of the lHodge &: KennicuttI (jl983l ) data set), S/^j is no longer detected in NGC 6946 
by our observations, and T^hi is at least as large as the value for the Milky Way. We conclude that 
given the similarity of Nhii /^pc^ and T,hi between these two galaxies, we would expect that T,h^ 
in NGC 6946 should be also be > 0.1 M© pc'^ at 12 kpc. 

Complicating these comparisons is the uncertainty in the NGC 6946 distance. To date, there 
are no Cepheid-based distance estimates to NGC 694 6, which has a low G alactic latitude. Recent 
estimates are (i n Mpc): 5.7ifc 0.7 ba sed on type II SN (jEastman et al.lll996l ). 6. 4 ±0.4 based on the 



brigh test stars (jSharina et al.l 119971 ). and 5.9 it 0.4 based on blue supergiants (jKarachentsev et al 



2OOOI ). As with any comparison based on absolute scales, our comparison changes if the the 6 Mpc 
distance estimate is incorrect by more than ±0.5 Mpc. 

The mapping of the CO disk presented here is sensitivity- limited; our 3a detection limit is 
2-3 M.Q pc~^. A combination of low T^x and a low filling factor for molecular clouds in the beam 
combine to create a detection limit rather than a real threshold. There is also the possibility that 
declining metallicity in the outer disk lowers the CO emissivity relative to j although this effect 
has yet to be seen. There are probably clouds here, but more widely dispersed than they are in the 
inner parts of the NGC 6946. 

But our suspicion that there is substantial molecular gas at large galactic radii is not based 
entirely on a comparison to the Milky Wa y. There a.re theoretical reasons to expect molecular 



gas at radii beyon d our detected CO edge (lElmegreen 



198' 



1 



Honma et al. 1995). According to Elmegreen Sz Parravanol ( 1994 ) molecular clouds are not formed 



1993 



Elmegreen &: Parravanol 1 19941 : 



Caveat: these are not azimuthal averages, but rather represent volumes at the sun and a variety of viewing angles 
centered on a Galactic anticenter vector, and may not be representative of the Galaxy as a whole. 
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where the ISM pressure, Pism, becomes too low to facihtate the formation of cold, dense clouds. Or 
conversely, molecular clouds form in rGgions where Pism 

> Pmin- The outer HI disk of NGC 6946 
has considerable surface density structure that exists on kpc in-plane size scales and where Pism 
may locally exceed this minimum. To calculate the mid-plane gas pressure we use: 



P. 



0.5 vr G ^gas^ 



where we have ignored the stellar co ntribution to the pressure which can only act to increase the 
value ( Elmegreen &: Parravandll994l ). Along the outer (R~ 6', 10 kpc) HI gas arms, Pism 

~ 2 X 10^ 



cm~^ K which is close to Pism values at the edge of the CO disk. The outer gas arms are where one 
would expect to find molec ular gas, if we had deeper maps. Several of the HII regions identified by 
Hodge &: KennicuttI (|l983l ) are located beyond the edge of our CO disk, particularly to the south, 
and it is likely that they ar e associated with mol ecular clouds that are too beam-diluted and cold 
to show up in our images. iFerguson et al.l (jl998l ) have detected Ha emission associated with the 
formation of massive stars, up to 3' beyond our CO disk edge and precisely at the local maxima in 
the outer HI arms of NGC 6946. These star forming outer arms would be good places to look for 
CO with a smaller beam and higher sensitivity. 

Based on theoretical arguments and the detection of Ha in NGC 6946 at large galactic radii, 
the answer is: no, the "edge" of the CO disk at a galactocentric radius of 5' (~ 9-10 kpc) in our 
images is not a real edge, it is a detection limit. There is no significant change in properties to 
indicate that the disk is undergoing a change at this radius. There could easily be molecular clouds 
that escape detection beyond this radius 



6. Conclusions 

We have obtained deep CO(l-O) and C0(2-l) images of a 10' x 10' region centered on the Sc galaxy, 
NGC 6946, with 55" and 27" resolution using the NRAO 12 Meter Telescope. We have combined 
these deep CO images with VLA HI images to obtain images of the neutral gas in this galaxy. To 
summarize our findings: 

1) The CO(l-O) and C0(2-l) disks shows many of the same features seen in the optical disk: 
a strong nuclear peak of emission, molecular gas arms spatially coincident with the asymmetric 
optical arm pattern, and a 10' diameter inner disk filled with CO emission. The CO disk is roughly 
exponential with a scale length of 1.2' it 0.2' or 2 it 0.3 (D/6 Mpc) kpc. 

2) We obtain 3.3 x 10^ (D/6 Mpc)^ M© for the molecular hydrogen gas mass of NGC 6946. This can 
be compared to the total atomic hydrogen gas mass of Mhj = 7.0 x 10^ (D/6 Mpc)^ Mq (which 
includes a VLA missing flux estimate). Molecular hydrogen constitutes 1/3 of the interstellar 
hydrogen gas mass and 2% of the dynamical mass of NGC 6946. At the nucleus, Sj/j = 170 
Mq pc~^, the molecular gas surface density is > 20 Mq pc~^ along the strong optical arm pattern, 
and ~ 10 Mq pc~^ regions beyond. 

3) The mean value for ri2 = Ico{2-i) / Ico{i-o) the the nucleus and optical arm pattern is 
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~ 0.83 lb 0.14. This is consistent with optically thick gas, with an excitation temperature of 10-15 
K. Interarm CO shows the largest variation in ri2, with a range of 0.35 to 2. Only 4% of the CO 
ratio map has ri2 < 0.6 which is likely to be cold or subthermally excited CO (T^x < 4 K). Fully 
half of the CO disk of NGC 6946 has ri2 > 1- The high ri2 appears to indicate the presence of 
widespread, optically-thin gas in between the spiral arms. 

4) Molecular gas dominates the total gas surface density within a galactocentric radius of 3.5' (6 
kpc). CO emission fills the central depression in the HI disk. Beyond the nuclear region, there is 
a good correlation between the CO and HI gas, which both peak on the spiral arms. The inner 
disk arm pattern transitions smoothly from predominantly CO arms to predominantly HI arms at 
the outskirts of the CO disk. At the nucleus, T,gas = 240 Mq pc~^, the total gas surface density 
is > 20 — 45 M0 pc~^ along the optical arm pattern, and ^ 5 — 10 Mq pc~^ along the outer gas 
arms {'^gas = l-36(Sij2 -|- Hhi))- C0(2-1) emission peaks are coincident with clumps of HI, Ha, 
and FIR emission outside the nuclear region, presumably an indication of warmer molecular gas 
and dissociated H2. 

5) We find strong correlations between star formation tracers (except for Ha) and CO emission in 
NGC 6946. The correlation coefficients are ~ 0.95. When T,gas is used instead of just the molecular 
gas tracer in a radial comparison, the dependence on galactic radius is reduced. Star formation is 
more closely related to the total gas surface density than the molecular gas surface density alone, a 
Schmidt law relationship. This also implies star formation is likely well beyond the radial limits of 
our observed CO disk in regions of predominantly HI gas. Surprisingly, Ha is not as well correlated 
with CO or Yjgas as are the other star formation tracers. The correlation coefficient is significantly 
less, ~ 0.65, and the ratio formed from Ha and CO or Ugas declines more rapidly with galactic 
radius. 

6) The star formation efficiency in the disk of NGC 6946 is relatively uniform; 2.8 Gyr when 

expressed in terms of the gas consumption time-scale. This uniformity implies 'Ssfb ^ ^gas- The 
star formation rate at the nucleus is due to the high molecular surface density rather than an 
elevated star formation efficiency. 

7) The edge of the observed CO disk represents a detection limit rather than a threshold for 
molecular cloud formation. We suspect the exponential fall of the CO disk continues smoothly past 
our detection limit. The mid-plane gas pressure in the outer HI arm structure, Pism ~ 2 x 10^ 
cm~^ K is close to the value at the limits of the observed CO disk. These gas arms should support 
the formation of molecular clouds which could be detectable in higher resolution observations. 

The authors would like to thank the NRAO 12m telescope operators, Duane Clark, Paul Hart, 
Victor Gasho and Harry Stahl, for their help and company during the many observing sessions 
in which the data for NGC 6946 and other galaxies were acquired. We thank L. Tacconi giving 
for us permission to retrieve and use the HI 21 cm line observations from the VLA archive. And 
we appreciate the careful reading and critical suggestions made by the anonymous referee. This 
work was supported in part by NSF grants ASTOO-71276 and AST03-07950. We made use of the 
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Fig. 1. — CO(l-O) channel maps for NGC 6946. Grey scale ranges from to 0.54 K. Contours are 
at -0.070, 0.070, 0.14 and 0.21 K. The rnis noise in a channel with no emission is 0.033 K. The 55" 
(FWHM) circular beam is displayed in the -61 km s~^ channel at the lower right. Each channel is 
labeled with the channel velocity (LSR). The channels are separated by 10.4 km s~^. 



-25- 




Fig. 2. — CO(l-O) integrated and peak intensity maps for NGC 6946. a) CO(l-O) integrated 
intensity map. Grey scale ranges from to 50 K km s^^. The contour levels arc 1, 2.5, 5, 7.5, 
10, 12.5, 15, 20, 30, 50, and 70 K km s^^. The 55" (FWHM) circular beam is displayed at the 
lower left, b) CO(l-O) peak intensity map. Grey scale ranges from to 0.52 K. The contour levels 
are 0.10, 0.17, 0.23, 0.30, 0.37 and 0.43 K. c) DSS uncalibrated optical image with the same Iio 
contours shown in (a). 
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Fig. 3. — C0(2-l) channel maps for NGC 6946. Grey scale ranges from to 0.85 K. Contours are 
at -0.17, 0.17, 0.34 and 0.51 K. The rms noise in a channel with no emission is 0.085 K. The 27" 
(FWHM) circular beam is displayed in the -60 km s~^ channel at the lower right. Each channel is 
labeled with the channel velocity (LSR). The channels are separated by 10.4 km s~^. 
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Fig. 4. — C0(2-l) integrated and peak intensity maps for NGC 6946. a) C0(2-l) integrated 
intensity map. Grey scale ranges from to 50 K km s^^. The contour levels are 5, 10, 15, 30, 50, 
75, and 100 K km s^^. The 27" (FWHM) circular beam is displayed at the lower left, b) C0(2-l) 
peak intensity map. Grey scale ranges from to 0.8 K. The contour levels are 0.26, 0.43, 0.60 and 
0.77 K. c) DSS uncalibrated optical image with the same I21 contours shown in (a). 
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Fig. 5. — Azimuthal plots of CO in NGC 6946 at selected radii. The Ico maps were transformed 
to face-on images and the azimuthal variation in Iqo was measured at constant galactic radii. Due 
south is at ^ = 0°, west is at = 90°. C0(2-l) is traced by the dashed line, CO(l-O) by the solid 
line. The data are presented at their native beam sizes, 27" and 55", for the C0(2-l) and CO(l-O) 
respectively, a) The azimuthal variation at R = 2.5' (4.4 kpc). b) The azimuthal variation at R = 
3.5' (6.1 kpc). 
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Fig. 6.— The CO(2-1)/CO(1-0) ratio in NGC 6946. a) ria = hi/ ho- The contours are at ri2= 
0.8, 1.0 and 1.2. For the colored regions: red represents ri2 < 0.8; orange for 0.8 > ri2 > 1.0; 
yellow for 1.0 > ri2 > 1.2; green for ri2 > 1.2; blue for ri2 ~ 2. b) Digital Sky Survey optical 
image with ri2 in color. The color scheme roughly follows that of figure a) and is intended to show 
the relationship between the regions of high/low ri2 and the optical stellar disk. 
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Fig. 7. — CO and HI in NGC 6946. a) Ihi (grey) and /lo (contours). Grey scale ranges from 0.5 
to 3.7 Jy beam^^ km s^^. The Iio contours and are at 1, 2.5, 5, 10, 15 and 35 K km s~^. b) Ifji 
(grey) and I21 (contours). Grey scale ranges from to 3.7 Jy beam~^ km s~^. The I21 contours 
and are at 5, 10, 15, 30 and 50 K km s~^. 



-31 - 



2.5 



o 1-5 

Oh 



0.5 



1 1 1 1 1 

"\ \ 


1 1 1 1 1 1 1 1 




III 


\ \ 

- \\ 

- \\ 


^HI 




: \\ 
- \\ 


Ej^j + 1.2 




\ 

; \^ \ 


1.36 (S^, ^ 










1 1 1 1 1 


\ 

\ 

\ 

1 1 1 . 1 1 1 1 1 1 1 


III 



4 6 8 

R (arcmin) 



10 



12 



Fig. 8. — Gas surface density vs radius in NGC 6946. Plotted are: Sijj derived using the standard 
conversion factor; Ylui from the observed emission; Ylm increased for the VLA missing flux esti- 
mate; and the total gas surface density increased by a factor of 1.36 to include the He and heavier 
element content. All have been corrected for inclination, (cos i). At our assumed distance of 6 
Mpc, 1' = 1.75 kpc. 
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Fig. 9. — The neutral gas surface density in NGC 6946. a) A map of the total neutral gas surface 
density, corrected for inclination (cos(40°)) and increased by a factor of 1.36 for the He and heavier 
element content. This map does not include a missing HI flux estimate. Contours are at 1, 2, 4, 
6, 10, 20, 40, 60 and 100 M0 pc^^. b) A false color image with CO(l-O) in green and HI in red. 
Regions of nearly equivalent column densities show up in orange. The C0(2-l) emission has been 
added in blue to bring out the bright CO nucleus. 
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Fig. 10. — CO and star formation tracers in NGC 6946. The CO(l-O) and C0(2-l) beam sizes are 



55^' an d 27" respectively, a) I21 in contours and Ha in grey. The Ha image is from [Ferguson et al 



(jl998l ) and is not presented here in caUbrated units. The Ha image was regridded to 2" pixels. 
The I21 contours are at 5, 10, 15, 30, and 50 K km s~^. b) I21 in grey with IRAS 60 //m contours. 
The grey scale ranges from to 40 K km s~^. The 60 /xm contours are at 10, 20, 30, 40, 60, 90, 
120 and 200 MJy sr'^. The 60 /im beam size (FWHM) is 45" x 41", pa = 21° and the rms noise 
level is 0.6 MJy sr^^. c) Iiq in grey with IRAS 100 ^m contours. The grey scale ranges from to 
40 K km s-^. The 100 //m contours are at 20, 40, 60, 80, 120, 160, 200 and 300 MJy sr'^. The 
100 ^m beam size (FWHM) is 69" x 65", pa = 20° and the rms noise level is 2 MJy sr~^. d) Iio 
in grey with 21 cm continuum contours. The grey scale ranges from to 40 K km s~^. The 21 cm 
contours are at 5, 7, 10, 12, 18, 22, 25, 36 and 48 mJy beam^^ The beam size size (FWHM) is 
49" X 42", pa = 73° and the rms noise level is 1.2 mJy beam~^. 
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Fig. 11. — Radial star formation rates and efficiencies in NGC 6946. a) '^sfr derived from the 
I RAS luminositi es, Lfir] non-thermal radio continuum lur ninosity, L ntr, using the SFR calibration 
of lCondonI (Il992l') . andLntc using the SFR calibration from lBehl ([2003); and as predicted by Schmidt 



law of 



KennicuttI (|l998al ) from T,gas- b) The star formation efficiency expressed in terms of tsfe. 



the time it would take to consume all of the molecular gas at the given star formation rate. The 
legends in both figures indicate the source of the SFR calibration. At our assumed distance of 6 
Mpc, 1' = 1.75 kpc. 
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Fig. 12. — Ratio of total gas surface density to 21 cm continuum in NGC 6946. The gas surface 
density maps was convolved to a 70" FWHM beam size and clipped at a 3a level prior to forming 
the ratios. The resulting ratio maps were sampled on a 35" grid. At our assumed distance of 6 
Mpc, 1' = 1.75 kpc. 
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Table 1. Global Properties of NGC 6946 



Property Value 

Hubble type ^ Scd 

R.A. (J2000) ^ 20'^ 34™ 51.^6 (±0.^4) 

Dec. (J2000) ^ 60° 9' 8.7" (±1.6") 

l^\h^^ 95.7°, 11.7° 

Distance ^ 6 Mpc 

v^s^ 50.5 (±1.3) km s'^ 

Inclination'^ 40° (±10°) 

Position Angle 242° (± 1° ) 

Vraax 170 (±40) km S'^ 

RVrr^a. " 5.l'(±0.7') (9.0 kpc) 



vmax 

d 



n 1.07 (±0.07) 

l0 



MH2totar 3.3x109Mo 



Mm observed ^ 5 X 10^ M© 

Mm ^ 7 X 10^ M© 

Mgas ^ 1.4 X IQio M0 

Mdyn ' 1.9 X lOii Mo 



^obtained from the NASA/IPAC Extragalactic Database. 

^dynamical center based on fit of Brandt rotation model to the HI 
data. 

'^Sharina et al. (1997). 

'^kinematic parameters from fit of Brandt rotation model to the 
HI data. 

^detected, using a 2.0 x 10^° K km s~^conversion factor, 
detected, no VLA missing flux estimate, 
^includes VLA missing flux estimate. 

'^includes factor for He and heavier elements: 1.ZG{Mh2 + Mhi) ■ 
Hotal dynamical mass based on Brandt model fit parameters. 
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Table 2. Properties of 3 Late-type Galaxies and the Milky Way 



Cxfil a.xv 


IC342 ^ 


NGC 6946 M83 ^ 


Milkv Wav 


Hubble type 


Scd 


Scd 


SBc 


SBbc 


Distance (Mpc) 


3.3 


6 


4 




(arcminutes) ^ 


20 


11 


11 






0.23 


0.38 


0.15 






1.1 


1.0 


1.0 




^dynamical (W^^ Mq) 


1.8 f 


1.9 


0.7 


~ 3 e 


MH.total (10^ M0) 


16 


33 


23 ' 


25 j 


Mmtotai (10« Mo) k 


57 


70 


62 


48 1 


MhJMhi 


0.27 


0.47 


0.37 


0.52 


Mgas/Mdyn ™ 


0.055 


0.074 


0.16 


0.03 


nucleus ^y^Q pc'^) ° 


140 


170 


220 


400 ° 




11 


13 


23 


3.4 q 



^data from Crosthwaite et al. (2000,2001). 
'^data from Crosthwaite et al. (2002). 

'^latest distance estimate based on Cepheid variables Saha et al. (2002). 
'^using the radius at which the azimuthal average S_ff2 ^ 0.5 pc~^. 
^D25 from RC2, de Vaucouleurs, de Vaucouleurs & Corwin (1976). 

r 

dynamical mass from Brandt model using 3.3 Mpc adjusted from 2 Mpc. 
S within 20 kpc, Nikiforov et al. (2000), from rotation curve modeling, 
'^adjusted to 3.3 Mpc distance from 2 Mpc. 
'rescaled to a 2.0 x 10^° K km conversion factor. 
JCombes (1991). 

'^includes VLA missing flux estimates. 
'Kulkarni & Heiles (1987). 

°^gas includes factor of 1.36 for He and heavier elements. 
"^Mq pc~^ in 55' beam centered on the nucleus. 
°R < 400 pc, Scoville & Sanders (1987). 
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Pmean Ej^j from observed CO disk excluding the nucleus. 
I4OO pc < R < 14 kpc, Scoville & Sanders (1987). 
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Table 3. Radial Slope of Gas to Star Formation Tracer Ratios 



Gas Tracer ^ 


SF Tracer ^ 


Slope (Uncertainty) ^ 


-'21 


inn iim 
iuu yum 


n 1 Qfi ^'n ni 


hi 


60 /xm 


-0.189 (0.015) 


hi 


21 cm continuum 


-0.170 (0.014) 


hi 


Ha 


-0.224 (0.015) 


ho 


100 fim 


-0.102 (0.008) 


ho 
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^Maps wore convolved to a 70" circular beam size and clipped 
at 3a. The region containing the bright 21 cm continuum emission, 
tentatively identified as a background galaxy, at a = 20''34'"24.9^ 
5 = 60° 10 38 (J2000) was clipped from all the comparison maps. 

^Slope (ratio/arcminute) and rms uncertainty from a least squares 
fit to the ratio of gas tracer to formation tracer as a function of 
galactic radius, log{Gastracer /SFtracer) = A ■ r + B. The maps were 
sampled on a grid of half beam width points (35"). 

''The molecular gas surface density derived from CO plus the 
atomic gas surface density derived from HI, increased by a factor 
of 1.36 to account for He and heavier elements. 



-40- 



Table 4. Comparison of the Milky Way and NGC 6946 at Large R 





Milky ^'^ 
Way 1 


B kpc 


12 kpc 


NGC " 

6946 8 kpc 12 kpc 


Af/////kpc2 




21 


> 1.4 


~ 5 d ^ 1 d 


EhMq PC"') 




2.9 


> 0.15 


2.3 ± 1.7 .. ^ 






2.9 
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^using R© = 8 kpc. 

''values from discussion in Wouterloot et al. (1988) 
'^using distance to NGC 6946 = 6 Mpc. 

'^the HII region positions from Hodge & Kennicutt (1983) were sampled on a grid of 

points separated by 10". The points were regrouped into 20" wide galactic radius bins. 
The 12 kpc value was extrapolated from the resulting Nhii /^pc^ curve. 

^beyond our measured CO disk. 

r 

including the estimated flux missing from the VLA map would increase these by 1.2 
Mo pc-^ 



